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Lp-discrepancy of the symmetrized van der Corput 

sequence 
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Abstract 

It is well known that the Lp-discrepancy for p E [l,oo] of the van der Corput 
sequence is of exact order of magnitude 0{(logN)/N). This however is for p E 
(l,oo) not best possible with respect to the lower bounds according to Roth and 
Proinov. For the case p = 2 it is well known that the symmetrization trick 
due to Davenport leads to the optimal -L 2 -discrepancy rate 0{^/\og^V/N) for the 
symmetrized van der Corput sequence. In this note we show that this result holds 
for all p E (1, oo). The proof is based on an estimate of the Haar coefficients of the 
corresponding local discrepancy and on the use of the Littlewood-Paley inequality. 

Keywords: Lp-discrepancy, van der Corput sequence, Davenport’s reflec¬ 
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MSC 2000: 11K38, 11K31 

1 Introduction and Statement of the Result 

For an infinite sequence S = {xn)n>o of points in [0,1) the local discrepancy of its first 
N elements is defined as 

^ N-l 

DN{S,t) = JjY. -h 

where, throughout this paper, 1/(0;) denotes the indicator function of the interval I C 
[0,1]. The Lp-discrepancy for p E [1, oo] of S is defined as the Lp-norm of the local 
discrepancy, thus, for p E (1, oo) 

Lp,^(5) := ||n^,(5,-)||L, = (^£\D^{S,t)\Pdty\ 

For p = oo we have 


Loc,,n{<S) := ||n7v(*S, •)|]z,^ = sup |n7v(*s,t)|. 

tG[0,l] 

The Lp-discrepancy is a quantitative measure for the irregularity of distribution of a 
sequence modulo one, see, e.g., P, [T3] . It is also related to the worst-case integration 
error of a quasi-Monte Carlo rule, see, e.g., pm- 
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It is well known that for every p G (1, cxo] there exists a positive nnmber Cp with the 
property that for every sequence S in [0,1) we have 

log N 

Loo,n{<S) > Coo jy for inhnitely many N eN (1) 

and, for p G (1, cxo), 

\/log N 

Lp,Ni'S) > Cp ——— for inhnitely many N eN, (2) 

where log denotes the natural logarithm and where N denotes the set of positive in¬ 
tegers {1, 2, 3,..The result for p = cx) was shown by Schmidt [21] and the result for 

p G (1, oo) was shown by Proinov [20] based on results of Roth [23] and Schmidt [23] . 

Both lower bounds ([T]) and ([2]) are optimal in the order of magnitude in N. 


A prototype of a sequence with low discrepancy is the van der Corput sequence (in 
base 2). Let <p(n) denote the radical inverse of n G No in base 2 (where No = N U {0}) 
which is dehned as <p(n) := ^.,2“*“^ whenever n has binary expansion n = Yl!i=o nj2*, 

where n* G {0,1} for all i G {0,..., k}. The van der Corput sequence (in base 2) is the 
sequence V = {yn)n>o where Un = <p(n) for n G Nq. 

For the van der Corput sequence it is known that (see, e.g, BE]) 


lim sup 

N^CO 


NLoo^n{V) 

logiV 


1 

3 log 2 


and hence Too, 7 v(V) is of order of magnitude 0{(\ogN)/N) which is best possible in N 
according to ([T]). However, for p G (1, oo) it is known that (see, e.g., [H [22] for p = 2 
and [12] for general p) 


lim sup 


NLp^NjV) 

\ogN 


1 

6 log 2 ’ 


This means that Lp^jqiy) for p G (1, oo) is only of order of magnitude 0{{\ogN)/N) 
which is not best possible in N if we compare with (121) . 

One way out to overcome this defect of the van der Corput sequence is based on 
symmetrization which was initially introduced by Davenport for (nQ!)-sequences (see [H 
Theorem 1.75]). This method is also known as Davenport’s reflection principle. 

We dehne the symmetrized van der Corput sequence (in base 2) = {zn)n>o as 




ip[m) if n = 2m, 

1 — <p(m) if n = 2m -|- 1. 


Then it is known, see e.g. [31[H1[I2112I], that the L 2 -discrepancy of the symmetrized 
van der Corput sequence is of optimal order of magnitude in N compared to the lower 
bound in ([2]), i.e.. 


T2,7v(V^y“) « 


V^ogN 

7^ 


(3) 


Here and throughout the paper, for functions /, p : N —)■ M’*', we write g{N) -C f{N), if 
there exists a C > 0 such that g{N) < Cf{N) for all G N, > 2. If we would like to 
stress that C depends on some parameter, say p, this will be indicated by writing -Cp. 

It is the aim of this paper to show that the estimate ([3]) holds for all p G (1, cxo). We 
show: 
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Theorem 1 For every p G (1, cxo) we have 




y/hgW 

7A~ 


The proof of this result is based on the Haar function system (in base 2) and will be 
given in Section [21 First we collect some auxiliary results in the following section. 

2 Auxiliary Results 

In order to estimate the Lp-discrepancy of we use the one-dimensional Haar system. 
Haar functions are a useful and often applied tool in discrepancy theory, see e.g. mm 

mmm- 

To begin with, a dyadic interval of length 2“-^, j G Nq, in [0,1) is an interval of the 
form 



for m = 0,1,..., 2-^ — 1. 


We also dehne /_i,o = [0,1). The left and right half of / = Ij^m are the dyadic intervals 
/■'“ = Am = Ij+i, 2 m and I~ = respectively. The Haar function 

hi = hj^rn with support I is the function on [0,1) which is -|-1 on the left half of /, —1 
on the right half of I and 0 outside of I. The Loo-normalized Haar system consists of 
all Haar functions hj^m with j G No and m = 0,1,..., 2-^ — 1 together with the indicator 
function h_i_o of [0, !)• Normalized in T2([0,1)) we obtain the orthonormal Haar basis 


of T 2 ([ 0 , 1 )). 

The Haar coefficients of a function / G Tp([0,1)) are dehned as 



where here and later on we use the abbreviations N_i := NoU{ —1}, := {0,1,..., 2^ — 

1} for j G No and D_i := {0}. 

In the following, we will compute the Haar coefficients of the local discrepancy of 


i.e.. 



Proceeding the computation of the Haar coefficients, we collect some properties of the 
radical inverse function (p{n) which we will need in the proof of the essential Lemma O 

Lemma 1 The following relations hold for the radical inverse function p: 


1. ‘pAs) = ^p{s) for allj,s G Nq, 

2. pAtA)) — P A 3 ^ CLf^d m G {0,..., 2-^ — 1}, 

3. pin) G Ij^m if o,nd only if n = 2ipim) + 2-^s for some s G No, 

I i< Elo |1 - ^Tis)\ <j + lforallAe No- 

Proof. 1. Let s = Ei=o -5*2*, where s* G {0,1} for alH G {0,..., k}. Then p{2is) = 

TiJ:loS^2^A = Eto = 2-^ EtoS*2—1 = X^is). 
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2. Since 0 < m < 2-^ — 1, m has a binary representation of the form m = I]i=o 

where rrij G {0,1} for alH G {0, ..., j — 1}. Then 2^if{m) = I]i=o and 

therefore = Eto = 2-^ Eto "^*2* = m2-K 

3. We write n in the form n = n + 2^s, where h G {0,1,..., 2-^ — 1} and s G Nq- Then 

(p{n) = n'2~^ for some n' G {0,1,..., 2-^ — 1} as one can verify easily. We have 
(p(n) = ^ + We see that ^p{n) G Ij^m is true if and only if n' = m. But from 
Point [21 we know ip~^{m2~^) = 2^ip{m) and the proof of Point [31 is done. 


4. To begin with, we verify the relation |1 — 2(p(2n)| + |1 — 2ip{2n + 1)| = 1 for all 
n G Nq. Therefore we observe that (p(2n) < | for all n G Nq. We also have 
(p{2n + 1) = (p{2n) + (p(l) = (p(2n) + hence 


1 - 2ip{2n)\ + |1 - 2ip{2n + l)\ 


=1 — 2ip{2n) + 


1 - 2 {ip(2n) + 


=1 — 2Lp{2n) + I — 2Lp(2n)\ = 1. 


This leads to 

A {A-1)/2 (A-1)/2 . 

^|l-2(p(s)|= ^ {|l-2(p(2n)| + |l-2(p(2n + l)|} = ^ 1 = 

s=0 n=0 n=0 ^ 

for odd A and 

A (A-2)/2 

^|l-2^(s)|= 5: {|l-2^(2n)| + |l-2^(2n + l)|} + |l-2v9(24)| 

s=0 n=0 

(A-2)/2 ^ 

= ^ 1 + |1 - 2(p(y4)| = — + |1 - 293(y4)| 

n=0 ^ 

for even A. 

Hence in both cases we have ^ < E^=o |1 ~ 2(p(s)| < ^ + 1. 

□ 


Lemma 2 The Haar coefficient /i-i,o of the local discrepancy ■) satisfies 


1 /^- 1,0 



‘fijAt) ^ J_ 
N — 2N 


if N = 2M, 
ifN = 2M + l. 


Proof. We have 




7V-1 




n=0 ' 
7V-1 


'i-[0,t){Xn) dt 


N-1 


jY X/ o o /\r ^ 


n=0 


N 


Xr. 


n=0 


1 

2 
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We therefore have to investigate the sum Y^n=o ^n- H N = 2M, then we have 

2M-1 M-1 M-1 jY 

Xn= ^ (1 - ^{m)) = M =—. 

n=0 n=0 m=0 

For N = 2M + 1, we find 

2M 2M-1 N — 1 

+ ^{M) = M + ^p{M) = —--h ^ (M). 

n=0 n=0 ^ 

This leads to the desired result. □ 

Remark 1 Lemma [2] is a crucial fact. It shows how the symmetrization trick keeps 
the Haar coefficient /i-i,o small enough in order to achive the optimal Lp-discrepancy 
rate. Let us compare this with the behavior of the Haar coefficient of the local 

discrepancy of the first N terms of the usual (not symmetrized) va,n der Corput sequence: 
for <N < 2"*+^ it is known that 


m—1 

i+E 

r=0 

where ||x|| denotes the distance of x to the nearest integer, see, e.g., [3 Proposition 1]. 
Further we know from m Theorem 3] that 




m—1 


max 


E 


N 




Hence it follows that there exists some constant c > 0 such that 

P-i'fo ^ ^— for infinitely many N ^N. 

Therefore, for the usual van der Corput sequence, already the first (and only the first, 
c.f. Lemma E]) Haar coefficient of the local discrepancy has the “bad” order of magnitude 
(logiV)/iV. 


In the following, we write := f) and denote by D%{t) the local 

discrepancy of the first N elements of the sequence {(p{n))n>o and analogously D](~‘^{t) 
the local discrepancy of the first N elements of the sequence (1 — (fi(n))n>o- Let 
and be the Haar coefficients of these three functions. 

Lemma 3 For N = 2M we have 

D7^7m = I (out) + D'-^nt)) 


and for N = 2M + 1 

C?M+,(() = ((A^ + 1)0?,+,(i) + . 
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Proof. For N = 2M we have 


1 /M-l M-1 \ 

^2M (^) E 1[0,4)(¥^W) + E l[0,i)(l - Pi^)) - t 

\n=0 n=0 / 

1 / 1 AT-l M-l 

= 2 (^M ^ E l[0,h(l - V?W) - t 

=i {Dt,{t) + Dl-^(t )). 


and for N = 2M + 1 we obtain 

1 / M M-l \ 

^2M+l(^) E l[0,t)(v?(^)) + E l[0.i)(l - Pi.'^)) - t 

ZIVI + i j 

1 / M M-l 

= 9 aJTT E l[0,b(v^W) -{M + l)t+ Y, l[ 0 ,t)(l - p{n)) - Mt 

ZlVl + i 


□ 


Corollary 1 IFe have 


\P 


N,syra i 
j,m. \ 


< 


2aEt 


+ M\ia 


j,‘m 


I) 


if N = 2M, 
ifN = 2M + l. 


Proof. We consider linearity of integration and the triangle inequality to obtain the 
result from Lemma [21 □ 


We proceed with the calculation of laj^m in the case j G Nq and hrst prove the 
following general lemma. 

Lemma 4 Let j G Nq and m G ©j. Then for the volume part f{t) = t of the discrepancy 
function we have 

N.M) = 

and for the counting part g{t) = Y.n=o l[o,q(2^n) we have 

2^—j — l N—l 

Tj,m{9) = E (|2m +1 - - 1), 

n=0 

o 

j,m 

where Ij,m= (p, denotes the interior of Ij^m- 
Proof. Of course, 

rl rm2~^ p{m-\-l)2~^ 

Pj,m{f) = / t hj^rriit) dt = / tdt- fdf = -2“^^"^ 

Jo Jm 2 -i Jm 2 - 2 + 2 - 1-1 

The Haar coefficients of g are given by 

/•i / X \ 1 /■! 

Tj,m{g) = I E d[ 0 ,t){Xn)hj^m(t) j df = — ^ / l[o,q df. 

“'0 V^'' n=0 / n=0^ _^. 

In 
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We analyse X„. If Xn ^ Ij,m or Xn = it is evident that X„ = 0. One can check by 
simple integration, that in the case Xn G we have X„ = m2“t — Xn, and if Xn G 
then X„ = —— (m2“t — Xn). These resnlts can be combined to 

In = + 1 - 2^~^^Xn\ - 1) if Xn elj,m ■ 

The claimed resnlt follows. □ 


Now we are ready to show a central lemma. 
Lemma 5 We have 


I/it' 

for all 0 < j < |'log 2 N~\ and 

L.Nv’I ^ I N,l-ip, ^ n-2j-2 
\r^j,m \ \H'j,m \ ^ 

for all j > flog 2 N]. 

Proof. We start with Xn = piji) and investigate the snm 


1 1 


A^,l —(/J I 


I — 


< 


1 1 
N 2t 


N-l 


Y, (|2m +1 - 2^+VHI - 1), 


n=0 

ipin)eij 


which, according to Lemma [H we can transfer to 

A 

Y ( 2m + 1 — 2^~^^Lp (2^(p{m) + 2^s^ — 

S = 1 ^ 

= |;(| 2 m+l- 2 »+‘ Q + AVs) 


A A 

= i: (ii - 2'+v(2's)i -1 )=j: (ii - 2^(s)i -1). 

s=l s=0 

We still have to specify the npper index. The conditions 0 < n < iV — 1 and n = 
2t(p(m) +2ts lead to s < — (p(m), this is why we choose A := \/d^ — Lp{rn)\. We also 

O 

see that there are no elements of {xq, Xi, ..., xat-i} contained in if 2^ip{m)+2^ > N, 
which is fnlhlled if 2^ > N or j > [log 2 iV]. Regarding that fact, we immediately 
conclude from Lemma 0] that \fifA\ ~ 2 “^t -2 fQj, j > |■log 2 iV]. We proceed with the 
case j < [log 2 A^] and have 

2-t-i ^ 


An = -1^ E (|1 - 2v>(i>)l - 1) + 2-"'-" 


< 


N s=o 
2-t-i (A 


N 
2-i-2 

N 


+ 1 - (A + 1)) + 2-2t-2 


A + < 


2-J-2 /N-l 


N 


2t 


- </?(m) - 1) + 2 


>-2j-2 
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We also find 


i (2-«-2 + (^(m) + < 1 (2-2J-2 + 2-i-') . 


A.':r>^(^-(-4+i))+2-«-^ 


N 

2-^-1 //I 


AT 


\ 9 , 2 ^ / iV — 1 (jjfm) \ n_- n 

+ 1 + >-+ 1 + 

2 / - iV V 2^+1 2 / 


= ^ 2 


-2i-2 


+ 


' ip{m) 


-12 


-3- 


N \ ' V 2 

By combining these results we finally obtain 


1 > A (2-^j-^ - 2-^-1 


\l4^\ < ^ (2-=^-= + 2-J-‘) < ^ ^ 


N 21 


as claimed. 


We turn to the estimation of l/i^^ which can be treated similarly to |/i 
begin with, we observe that 

N-l 

Y, {\2m + 1 - - ip{n))\- l) 


N,ip< 
■j,m I 


To 


n=0 

o 


N-l 


Y (|2m + l-2^+^(l-y?(n))|-l) 


n=0 
o 

B 


Y (|2m + 1 - 2^+^ (l - V? (2V(2^' - m - 1) + 2ls)) I - l) 


2m + 1 - 2^'+^ 1 


+ 4>{2^s) 


S = 1 

B 

= E 

s=\ 

B B 

= Y{\-^ + Ms)\-i) = Yi\^-Ms)\-i) 


2l - m - 1 
21 


S=1 


s=0 


In this expression, B := — (p{2^ — m — 1)J, which we deduce in the same way as 

the upper index A above. Completely analogously as above, we obtain 

I wp-yji ^ J_ (^2-21-2 I 2“1“1^ < — — 

for i < [log 2 A^]. The case j > [log 2 A^] also follows the same lines as above. □ 


Corollary 2 The Haar coefficients of the symmetrized van der Corput sequence for 
j e No fulfil 


Ih 


A^,sym I 

‘j,m I 




= 2-2^-2 ifj> nOg2iV]. 
Proof. We combine Corollary [T] and Lemma E] to obtain the result. 


□ 
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3 The Proof of Theorem [T] 


We are ready to show that the Lp-discrepancy of the symmetrized van der Corput 
sequence has optimal order in N for any p G (l,cx)). We apply the Littlewood-Paley 
inequality which involves the square function 


s{f) ■■= 




mGOj 


1/2 


of a function / G -hp([0,1)). 


Lemma 6 (Littlewood-Paley inequality) Let p G (l,cxo). Then there exist Cp,Cp > 
0 such that for every f G Lp([0,1)) we have 

cp||/IU,<||E/)IU,<c'p||/|U,. 


Proofs of these inequalities and further details also yielding the right asymptotic 
behavior of the involved constants can be found in [21 |2H1 [22]. Littlewood-Paley theory 
has already been used in the context of discrepancy before, see, e.g., ^mmm- 
Now we can give the proof of Theorem [H 

Proof. Throughout this proof, we simply write instead of Using Lemma E] 

with / = •) we have 


LpAV^n = ll^iv(V^^“,-)||L, 

«p \\s{D^{v^y^,-mL, 


1/2 


E E 

//SN-1 m&Bj 

y- 22max{0,i} V 1 , 

/SN_i mSOj 


1/2 


^p/2 


< E 22 max{0,j} 

ii6N_i 


^ y Lj,m 
mSitt,' 


1/2 


"p/2/ 


where we used Minkowski’s inequality for the Lp/ 2 -norm. Hence, in order to prove the 
result it suffices to show that 


^ ^ 2^2ma.x{0J} 

jeN-i 


Ep 


3 ,rn 


mG'Dj 



logiV 


Now Lemma [2] and Lemma E] give 


E 2 


2 max{0,j} 


jGN-i 



1 

4 

[log 2 N]-l 



< 

“4iV2 

|l[0,l)IUp/2 + Jp 

E 

y 1 / 

m^'Dj 

Lp/2 

j=0 

m£Bj 
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— s 


2-2j 


E 

m£B)j 


-p/2 


1 4 1 1 logiV 

+ 1 ) + ^ 

where we regarded the fact that 1/j m = 1 for a hxed j G Nq. The proof is 

complete. □ 
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